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Abstract 

We argue that the observation of a sizable direct CP asymmetry A^ S7 in the 
inclusive decays B —* X s j would be a clean signal of New Physics. In the Stan- 
dard Model, A^ S7 can be calculated reliably and is found to be below 1% in 
magnitude. In extensions of the Standard Model with new CP-violating couplings, 
large CP asymmetries are possible without conflicting with the experimental value 
of the branching ratio for the decays B —* X s j. In particular, large asymme- 
tries arise naturally in models with enhanced chromo-magnetic dipole operators. 
Some generic examples of such models are explored and their implications for the 
semileptonic branching ratio and charm yield in B decays discussed. 
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1 Introduction 



Studies of rare decays of B mesons have the potential to uncover the origin of CP viola- 
tion, which may lie outside the Standard Model of strong and electroweak interactions. 
The measurements of several asymmetries will make it possible to test whether the CKM 
mechanism of CP violation is sufficient, or whether additional sources of CP violation are 
required to describe the data. In order to achieve this goal, it is necessary that the the- 
oretical calculations of CP-violating observables in terms of Standard Model parameters 
are, at least to a large extent, free of hadronic uncertainties. This can be achieved, for in- 
stance, by measuring time-dependent asymmetries in the decays of neutral B mesons into 
particular CP eigenstates. In many other cases, however, the theoretical predictions for 
direct CP violation in exclusive B decays are obscured by large strong-interaction effects 
|I|]-pJ, which can only partly be controlled using the approximate flavour symmetries of 
QCD[|. 

Inclusive decay rates of B mesons, on the other hand, can be reliably calculated 
in QCD using the operator product expansion. Up to small bound-state corrections 
these rates agree with the parton model predictions for the underlying decays of the b 
quark M— [EJ. The possibility of observing mixing-induced CP asymmetries in inclusive 



decays of neutral B mesons has been emphasized in Ref. ||10|| . The disadvantage that the 
inclusive sum over many final states partially dilutes the asymmetries is compensated 
by the fact that, because of the short- distance nature of inclusive processes, the strong 
phases are calculable using quark-hadron duality. The resulting CP asymmetries are 
proportional to the strong coupling constant ^(mj). The purpose of the present paper 
is to study direct CP violation in the rare radiative decays B — > X s 7, both in the 
Standard Model and beyond. These decays have already been observed experimentally, 
and copious data samples will be collected at the B factories. As long as the fine structure 
of the photon energy spectrum is not probed locally, the theoretical analysis relies only 
on the weak assumption of global quark-hadron duality (unlike the hadronic inclusive 
decays considered in Ref. |TI|). Also, the leading nonperturbative corrections have been 
studied in detail and are well understood [[Ll[]-[[l^ . 



We perform a model-independent analysis of CP- violating effects in B — > X s ^ decays 

'cfl 



in terms of the effective Wilson coefficients Cj = C^^rrib) and C$ = Cg ff (mft) multiplying 



7 = ^ BLO^bn , Os = y -f± n^GTbR (1) 



the (chromo-) magnetic dipole operators 

em b _ g s m b 

in the low-energy effective weak Hamiltonian []19| |. We will allow for generic New Physics 
contributions to the coefficients C 7 and C 8 , possibly containing new CP- violating cou- 
plings. Several extensions of the Standard Model in which new contributions to dipole 
operators arise have been explored, e.g., in Refs. pD[-P5[. We find that in the Stan- 



dard Model the direct CP asymmetry in the decays B — > X s 7 is very small (below 1% 
in magnitude) because of a combination of CKM and GIM suppression, both of which 
can be lifted in extensions of the Standard Model. If there are new contributions to 
the dipole operators with sizable weak phases, they can induce a CP asymmetry that is 
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more than an order of magnitude larger than in the Standard Model. We thus propose a 
measurement of the inclusive CP asymmetry in the decays B — > X s 7 as a clean and sen- 
sitive probe of New Physics. For simplicity, we shall not consider here the most general 
scenario of having other, non-standard operators in the effective Hamiltonian. However, 
we will discuss the important case of new dipole operators involving right-handed light- 
quark fields, which occur, for instance, in left-right symmetric models. The interference 
of these operators with those of the standard basis, which is necessary for CP violation, 
is strongly suppressed by a power of m s /mb] still, they can give sizable contributions to 
CP-averaged branching ratios for rare B decays. 

Studies of direct CP violation in the inclusive decays B — > X s 7 have been performed 
previously by several authors, both in the Standard Model |24| and in certain extensions 
of it |ES[ |26|| . In all cases, rather small asymmetries of order a few percent or less 
are obtained. Here, we generalize and extend these analyses in various ways. Besides 
including some contributions to the asymmetry neglected in previous works, we shall 
investigate in detail a class of New Physics models with enhanced chromo-magnetic 
dipole contributions, in which large CP asymmetries of order 10-50% are possible and 
even natural. We also perform a full next-to-leading order analysis of the CP-averaged 
B — > Xg'j branching ratio in order to derive constraints on the parameter space of the 
New Physics models considered here. For completeness, we note that CP violation has 
also been studied in the related decays B — > X s [^?J , which however have a much 



smaller branching ratio than the radiative decays considered here. 



2 Direct CP violation in radiative B decays 

The starting point in the calculation of the inclusive B — > X s 7 decay rate is provided by 
the effective weak Hamiltonian renormalized at the scale /i = m b ||19| . Direct CP violation 
in these decays may arise from the interference of non-trivial weak phases, contained in 
CKM matrix elements or in possible New Physics contributions to the Wilson coefficient 
functions, with strong phases provided by the imaginary parts of the matrix elements of 
the local operators of the effective Hamiltonian J28j. These imaginary parts first arise 
at 0(a s ) from loop diagrams containing charm quarks, light quarks or gluons. Using 
the formulae of Greub et al. for these contributions |29 |, we calculate at next-to-leading 



order the difference Ar = T(B — > X s ~f) — T(B — > X^) of the CP-conjugate, inclusive 
decay rates. The contributions to Ar from virtual corrections arise from interference of 
the one-loop diagrams with insertions of the operators 2 and Og shown in Figure [TJ(a) 
and (b) with the tree-level diagram for b — > 37 containing an insertion of the operator 
O7. Here 2 = s^^q^qi^^b^ with q = c,u are the usual current-current operators in 
the effective Hamiltonian. We find 

G 2 F m 5 b aa s (m b ) 

Z_il -i. 



virt 



187T 4 



|2 



J lm[v u v* t C 2 C;} - (J - z «(*)) lm[v c v*C 2 C* 7 ] - M- lm[C 8 C* 7 ] \ , (2) 
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(a) (b) (c) 

Figure 1: Diagrams for b — > s"f(g) yielding non-trivial strong phases that can contribute 
to the CP asymmetry. The crosses indicate other possible attachments of the photon. 
The numbers inside the squares indicate which operators are inserted. 



where v q = V* s V q b are products of CKM matrix elements, z = (nic/mb) 2 , and 

v(z) = (^ + \nz + \n 2 z- + (\n 2 z - y + - i Inz^z 2 + 0(z 3 ) . (3) 

There are also contributions to Ar from gluon bremsstrahlung diagrams with a charm- 
quark loop, shown in Figure [l](c). They can interfere with the tree-level diagrams for 
b — > s'-yg containing an insertion of O7 or 0$. Contrary to the virtual corrections, 
for which in the parton model the photon energy is fixed to its maximum value, the 
gluon bremsstrahlung diagrams lead to a non-trivial photon spectrum, and so the results 
depend on the experimental lower cutoff on the photon energy. We define a quantity S 
by the requirement that E 1 > (1 — 5)E™ a,x , i.e. 5 is the fraction of the spectrum above 
the cut. 1 We then obtain 

Ar brems = glr^f^ z b ^ 5) (i m [ Vc v* t c 2 c* 7 ) - i im{v cV ;c 2 c;]) , (4) 

where b(z, 5) = g(z, 1) — g(z, 1 — 5) with 



g(z, y) = 6(y - Az) { (y 2 - Ayz + 6z 2 ) ln(</f + JJ^l) - M \ 2g) Jl - - 



Az V Az J A V V 

(5) 

Combining the two contributions, dividing the result by the leading-order expression for 
(twice) the CP-averaged inclusive decay rate, 

F(B -> X sl ) + T(B -> X s - 7 ) = | Vt C 7 | 2 , (6) 

and using the unitarity relation v u + v c + v t = 0, we find for the CP asymmetry 



r(s^x s7 )+r(s^x s - 7 ) 



E J >(l-5)E™^ 



1 In the parton model £;™ ax = m&/2 depends on the quark mass and does not agree with the physical 
boundary of phase space. Later, we shall discuss how this problem is resolved by including the effects 
of Fermi motion. 
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MC2C7) - I [v(z) + b(z, 5)} Im[(l + e s )C 2 C* 7 ] 
- ~ lm[C 8 C* 7 ] + I b(z, 6) Im[(l + e s )C 2 C* 8 ] \ , (7) 



where 



^ = ^ = §^-A 2 (zr / -p)=O(10- 2 ). (8) 

V 7 ^ V t b 

In the last step, we have expressed e s in terms of the Wolfenstein parameters, with 
A = sin^c ~ 0.22 and p, r\ — 0(1). We stress that (|7|) is an exact next-to-leading order 
result. All numerical coefficients are independent of the renormalization scheme. For 
consistency, the ratios of Wilson coefficients Cj must be evaluated in leading-logarithmic 
order. Whereas the bremsstrahlung contributions as well as the C 2 -C$ interference term 
are new, an estimate of the C 2 -C 7 interference term has been obtained previously by 
Soares p4j , who neglects the contribution of the function b(z, 5) and uses an approxima- 
tion for v(z). The importance of the Cg-CV interference term for certain extensions of 



the Standard Model has been stressed by Wolfenstein and Wu p5| , and the first correct 



calculation of its coefficient can be found in Ref. [|26 . 

In the Standard Model, the Wilson coefficients take the real values C 2 ~ 1.11, 
C7 ~ —0.31 and Cs ~ —0.15. The imaginary part of the small quantity e s is thus 
the only source of CP violation. Note that all terms involving this quantity are GIM 
suppressed by a power of the small ratio z = (m c /mj,) 2 , reflecting the fact that there 
is no non-trivial weak phase difference in the limit where m c = m u = 0. Hence, the 
Standard Model prediction for the CP asymmetry is suppressed by three small fac- 
tors: a s (mb) arising from the strong phases, sin 2 #c reflecting the CKM suppression, 
and (m c /mb) 2 resulting from the GIM suppression. The numerical result for the CP 
asymmetry depends on the values of the strong coupling constant and the ratio of the 
heavy-quark pole masses. Throughout this work we shall take a s (mf,) ~ 0.214 (corre- 
sponding to a s (m z ) = 0.118 and two-loop evolution down to the scale m h = 4.8 GeV) 
and \fz = m c /mb = 0.29. The sensitivity of the next-to-leading order predictions for 
inclusive B decay rates to theoretical uncertainties in the values of the input parameters 
as well as to the choice of the renormalization scale and scheme have been investigated 
by several authors. Typically, the resulting uncertainties are of the order of 10%. Since 
a discussion of such effects is not the purpose of our study, we shall for simplicity assume 
fixed values of the input parameters as quoted above. With this choice we find 

4p; S mW « 1-54% [l + 0.15 6(2, 5)] 77, (9) 

where < b(z, 5) < 0.30 depending on the value of S. With rj m 0.2-0.4 as suggested by 
phenomenological analyses j30], we find a tiny asymmetry of about 0.5%, in agreement 



with the estimate obtained in Ref. |24[]. Expression (0) applies also to the decays B 
Xdj, the only difference being that in this case the quantity e s must be replaced with 
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the corresponding quantity 



V ud V ub 



p — %r\ 
— p + ir) 



0(1). 



(10) 



Therefore, in the Standard Model the CP asymmetry in B — > X^7 decays is larger by 
a factor — (A 2 [(l — p) 2 + rf})~ 1 « —20 than that in B -> X s 7 decays. Note, however, 
that experimentally it would be very difficult to distinguish between inclusive B — > X s 7 
and B — > X^ 7 decays. If only the sum is measured, the CP asymmetry vanishes (in the 
limit where m s = rrid = 0), since 

AT SM (B - X s7 ) + Ar SM (S ^ X d7 ) cx im[K^; b (v;* s ^ + ok)] = (11) 



by unitarity. This has also been pointed out in Ref. [f24 . 

One might wonder whether our short- distance calculation of the CP asymmetry in 
inclusive B — > X s ^ decays could be upset by large long-distance contributions to the 
decay amplitude mediated by the current-current transitions, which could spoil quark- 
hadron duality. The most important process is likely to be B — > X S V followed by virtual 
V — > 7 conversion, where V = J/ip for the h — ► ccs transition, and V = p°,u> for the 
b — > uus transition. Using vector-meson dominance to estimate these effects pi"| , we 
find that the largest contribution to the asymmetry is due to J ftp — > 7 conversion and is 
at most of order 1%, i.e. at the level of the prediction obtained using the short-distance 
expansion. Hence, we see no reason to question the applicability of the heavy-quark 
expansion to predict the inclusive CP asymmetry. 

From (0) it is apparent that two of the suppression factors operative in the Standard 
Model, z and A 2 , can be avoided in models where the effective Wilson coefficients C-j 
and Cg receive additional contributions involving non-trivial weak phases. Much larger 
CP asymmetries of 0(a s ) then become possible. In order to investigate such models, we 
may to good approximation neglect the small quantity e s and write 



\C 7 \ 



{a 27 (5)lm[C 2 C* 7 } + a 87 lm[C 8 C* 7 ]+a 28 (5)lm[C 2 C* s }}, (12) 



where 



a$(6) = a s (m 6 )|^ - y [v(z) + b(z,8) 



4? = -gOz s (m b ) 



27 



a s (m b ) zb(z,6) . 



(13) 



The superscripts indicate that these results are obtained in the parton model. The 
values of the coefficients shown in the left portion of Table [l] for three choices 

of the cutoff on the photon energy: 5 = 1 corresponding to the (unrealistic) case of a 
fully inclusive measurement, 5 = 0.3 corresponding to a restriction to the part of the 
spectrum above ~ 1.8 GeV, and 5 = 0.15 corresponding to a cutoff that removes almost 



5 



Table 1: Values of the coefficients aij (in %), without (left) and with (right) Fermi 
motion effects included 



5 


a (p) a (p) a (p) 

"27 u 87 a 28 

(parton model) 


a 27 a 87 a 28 £f n [GeV] 
(with Fermi motion) 


1.00 
0.30 
0.15 


1.06 -9.52 0.16 
1.17 -9.52 0.12 
1.31 -9.52 0.07 


1.06 -9.52 0.16 0.00 
1.23 -9.52 0.10 1.85 
1.40 -9.52 0.04 2.24 



all of the background from B decays into charmed hadrons. In practice, a restriction 
to the high-energy part of the photon spectrum is required for experimental reasons. 
Whereas the third term in flT2|) will generally be very small, the first two terms can give 
rise to sizable effects. Since afj has a rather weak dependence on 5 and 0$ has none, 
the result for the CP asymmetry is not very sensitive to the choice of the photon-energy 
cutoff. Assume now that there is a New Physics contribution to C-j of similar magnitude 
as the Standard Model contribution, so as not to spoil the prediction for the CP-averaged 
decay rate in (§), but with a non-trivial weak phase. Then the first term in (0) may 
give a contribution of up to about 5% in magnitude. Similarly, if there are New Physics 
contributions to CV and Cg such that the ratio C 8 /CV has a non-trivial weak phase, the 
second term may give a contribution of up to about 10% x \Cs/Cy\. In models with a 
strong enhancement of |C 8 1 with respect to its Standard Model value, there is thus the 
possibility of generating very large CP asymmetries in B — ► decays. The relevance 
of the second term for two-Higgs-doublet models, and for left-right symmetric extensions 
of the Standard Model, has been explored in Refs. []25], p6| . 

In our discussion so far we have neglected nonperturbative power corrections to the 
inclusive decay rates. Their impact on the rate ratio defining the CP asymmetry is ex- 
pected to be very small, since most of the corrections will cancel between the numerator 
and the denominator. Potentially the most important bound-state effect is the Fermi 
motion of the b quark inside the B meson, which determines the shape of the photon 
energy spectrum in the endpoint region. Technically, Fermi motion is included in the 
heavy-quark expansion by resumming an infinite set of leading-twist corrections into a 
nonperturbative "shape function" F(k + ), which governs the light-cone momentum dis- 



tribution of the heavy quark inside the meson [12|, [T3[ . The physical decay distributions 
are obtained from a convolution of parton model spectra with this function. In the pro- 
cess, phase-space boundaries defined by parton model kinematics are transformed into 
the proper physical boundaries defined by hadron kinematics. For the particular case 
of the coefficients ajj? (5) in (13|), where in the parton model the parameter 6 is defined 
such that Ery > ~(1 — 6)1711,, it can be shown that the physical coefficients ctij(S) with 
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E~ > 



8)rriB are given by [[33 



aij(5) 



m B -m b 

I dk + F(k, 

m B (l-S)-rn b 



,(P) 



m B (l-S) 
m b +k + 



m B -m b 

I dk + F(k, 

m B (l-S)-m b 



(14) 



This relation is such that there is no effect if either the parton model coefficient is 
independent of 5, or if the limit 5 = 1 is taken, i.e. the restriction on the photon energy 
is removed. Several ansatze for the shape function have been suggested in the literature 
For our purposes, it is sufficient to adopt the simple form 



12, 13 



F(k+) = N(1- x) a e {l+a)x ■ x=^-<l 



(15) 



where A = ms — m^. The normalization iV cancels in the ratio in (|14]). The parameter 
a can be related to the heavy-quark kinetic energy parameter // 2 = — Ai p4| , yielding 



/i 2 = 3A 2 / (1 + a). In the right portion of Table || we show the values of the coefficients 
ctij(5) corrected for Fermi motion, using the above ansatz with m& = 4.8 GeV and /i 2 = 
0.3 GeV 2 . We also give the physical values of the minimum photon energy, E™ m = 
~(1 — 8)rriB- The largest coefficient, a 87 , is not affected by Fermi motion, and the impact 
on the other two coefficients is rather mild. As a consequence, our predictions for the 
CP asymmetry are very much insensitive to bound-state effects, even if a restriction on 
the high-energy part of the photon spectrum is imposed. 



3 Next-to-leading order corrections to B — ► X s ^y 

In the next section we shall explore in detail the structure of New Physics models with 
a potentially large inclusive CP asymmetry. A non-trivial constraint on such models is 
that they must yield an acceptable result for the total, CP-averaged B — > X s ~f branching 
ratio, which has been measured experimentally. Taking a weighed average of the results 



reported by the CLEO and ALEPH Collaborations [|| |3| gives B(B -> X s7 ) = (2.5 ± 
0.6) x 10~ 4 . We stress that this value is extracted from a measurement of the high- 
energy part of the photon energy spectrum assuming that the shape of the spectrum 
is as predicted by the Standard Model. For instance, the CLEO Collaboration has 
measured the spectrum in the energy range between 2.2 and 2.7 GeV and applied a 
correction factor of 0.87 ± 0.06 in order to extrapolate to the total decay rate |37j (see 
Ref. p3| for a critical discussion of this treatment). 



The complete theoretical prediction for the B — > X s ^ decay rate at next-to-leading 
order has been presented for the first time by Chetyrkin et al. [Q. The result for the 
corresponding branching ratio is usually obtained by normalizing the radiative decay 
rate to the semileptonic decay rate of B mesons, thus eliminating the strong dependence 
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Table 2: Values of the coefficients ku (in %) with Fermi motion effects included 



5 E™ in [GeV] 


kn k<i2 k$$ &27 kis &28 &77 


0.90 0.26 
0.30 1.85 
0.15 2.24 


75.67 0.23 8.47 -14.77 9.45 -0.04 3.47 
68.13 0.11 0.53 -16.55 8.85 -0.01 3.86 
52.18 0.03 0.11 -13.54 6.66 +0.00 3.15 



on the fe-quark mass. We define 

T(B ^ X sl )\ E ^ >(1 _ s)E 



T(B -> X c eu) 



6a 



v t *v tb 



cb 



K NLO (S) 



(16) 



where f(z) 



1 



Iz + 8z a 



Viz 2 In z is a phase-space factor, and the quantity 



-^nlo(^) — | C*r 1 2 + 0(a s , l/m|) contains the corrections to the leading-order result. 
Using cT 1 = 137.036 g| and \V*V tb /V cb \ « 0.976 as in Ref. @, we get 



£ 7 >(1-<5)E" 



2.57 x 10~ 3 K NLO (5) x 



X c eu) 



10.5% 



;i7) 



From now on we shall assume the value B(B — > X c eu) = 10.5% for the semileptonic 
branching ratio and omit the last factor. The current experimental situation of measure- 



ments of this quantity and their theoretical interpretation are reviewed in Refs. j|0], PJ| . 
The general structure of the quantity -Knlo is 



K NLO (S)= £ k ij (5)Re[C i C*] + k^(6)Re[C { 7 L) C* 



(18) 



i,j'=2,7,8 



where kij(5) are known coefficient functions depending on the energy cutoff parameter 
5, and C7 is the next-to-leading order contribution to the Wilson coefficient C^(m b ). 
In the Standard Model w 0.48 ||. Explicit expressions for the functions kij(S), at 
next-to- leading order in a s and including power corrections of order 1/ml, can be found 
in Ref. |^3[, where we correct some mistakes in the formulae used by previous authors. 
(The corrected expressions will also be given in an erratum to Ref. p8| .) Contrary to 
the case of the CP asymmetry, the impact of Fermi motion on the partially integrated 
B — > X s 7 decay rate is an important one for values of S that are realistic for present- 
day experiments. In Table 0, we show the values of the coefficients fcy corrected for 
Fermi motion [j3~3"| , using again m b = 4.8 GeV and fi 2 = 0.3 GeV 2 for the parameters of 
the shape function. We quote the results for three choices of the cutoff on the photon 
energy: 5 = 0.9 corresponding to an almost fully inclusive measurement, and 5 = 0.3 
and 0.15 corresponding to a restriction to the high-energy part of the photon spectrum. 
The choice 5 = 1 must be avoided because of a weak, logarithmic soft-photon divergence 
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E™ m [GeV] Ef n [GeV] 

Figure 2: Theoretical predictions for the integrated B — > X s ^f branching ratio for 
various choices of the parameters mi, and u, 2 ; left: n, 2 = (0.30 ± 0.15) GeV 2 for fixed 
rrib; right: nib = (4.80 ± 0.15) GeV for fixed ratio fi 2 /A 2 . The data point shows the 
CLEO measurement. 



in the prediction for the total B — > X s 7 branching ratio caused by the term proportional 
to ks8{S). Note that with a realistic choice of the cutoff parameter 5 the coefficient kgs of 
the term proportional to | Cs | 2 in ( |18D becomes very small. This observation will become 
important later on. With our choice of parameters, we obtain in the Standard Model 
B(B — > X s ^j) = (3.3 ± 0.3) x 10" 4 for 5 = 0.9 [j33fl , in good agreement with the results 
obtained in previous analyses 0, |2|, 



In order to illustrate the sensitivity of our results to the parameters of the shape 
function, we show in Figure |2] the predictions for the Standard Model branching ratio as 
a function of the energy cutoff E™ m = ^(l—5)rriB- In the first plot, we keep rrib = 4.8 GeV 
fixed and compare the parton model result (gray curve) with the results corrected for 
Fermi motion, using /z 2 = 0.15 GeV 2 (short-dashed curve), 0.30 GeV 2 (solid curve), and 
0.45 GeV 2 (long-dashed curve). This figure illustrates how Fermi motion fills the gap 
between the parton model endpoint at m b /2 and the physical endpointQ at m B /2. In 
the second plot, we vary m;, = 4.65 GeV (long-dashed curve), 4.8 GeV (solid curve), 
and 4.95 GeV (short-dashed curve), adjusting the parameter /z 2 in such a way that the 
ratio /i 2 /A 2 remains fixed. For comparison, we show the data point B(B — > Xg'-f) = 
(2.04 ± 0.47) x 10~ 4 obtained by the CLEO Collaboration with a cutoff at 2.2 GeV |7 



The fact that in the CLEO analysis the cutoff is imposed on the photon energy in the 
laboratory frame rather than in the rest frame of the B meson is not very important 
for the partially integrated branching ratio |33| and will be neglected here. Obviously, 
there is a rather strong dependence of the partially integrated branching ratio on the 
value of the 6-quark mass. In particular, by choosing a low value of m;, it is possible 
to get agreement with the CLEO measurement without changing the prediction for 
the total branching ratio. The important lesson from this investigation is that the 



2 The true physical endpoint is actually located at [m 2 B — (ttik + m n ) 2 ]/2mB ~ 2.60 GeV, i.e. slightly 
below ms/2 i=s 2.64 GeV. Close to the endpoint, our theoretical prediction is "dual" to the true spectrum 
in an average sense. 
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theoretical uncertainty in the prediction for the integral over the high-energy part of the 
photon spectrum is significantly larger than the uncertainty in the prediction of the total 
branching ratio. So far, this fact has not been taken into account in the comparison of the 
extrapolated experimental numbers for the total branching ratio with theory. Ultimately, 
the theoretical errors may be reduced by tuning the parameters of the shape function to 
fit the measured energy spectrum; however, at present the experimental errors are too 
large to make such a fit meaningful |H| . Below, we shall perform our calculations for the 
case 5 = 0.3 corresponding to E™ m « 1.85 GeV, which is large enough to be realistic for 
near-future experiments, yet low enough to be sufficiently insensitive to the modeling of 
Fermi motion. As we have pointed out before, the results for the CP asymmetry depend 
very little on the choice of cutoff. 



4 CP asymmetry beyond the Standard Model 

In order to explore the implications of various New Physics scenarios for the CP asymme- 
try and branching ratio in B — > decays it is useful to express the Wilson coefficients 
C 7 = C 7 s (m b ) and C 8 = C 8 s (mb), which are defined at the scale mj, in terms of their 
values at the high scale mw Using the leading-order renormalization-group equations, 
one obtains 

16 8 1 16 8 

C 7 = r/i C 7 (m w ) + - (r?i - C s (m w ) + J2 h i V ai > 

C 8 = v%Cs(m w )+j2h i V ai , (19) 

8=1 

where r\ = a s (mw) / 'ot s {rrib) ~ 0.56, and hi, hj and a, are known numerical coefficients 



44, 45 . For the Wilson coefficients at the scale mw, we write 



C 7 (m w ) = ~ A(xt) + Cr(m w ) , 

C 8 (m w ) = - 1 - D(x t ) + C«™(m w ) , (20) 

where the first terms correspond to the leading-order Standard Model contributions 
4f|. The y 

are known functions of the mass ratio Xt = (m t (my/) / /r rnw) 2 \ which we 



evaluate with m t (m w ) w 178 GeV (corresponding to a pole mass of 175 GeV). This 
yields ^A(x t ) ~ 0.20 and ^D(x t ) ~ 0.10. Using a similar evolution equation for the 

next-to-leading coefficient C7 PS|, we find 3 



C 7 w -0.31 + 0.67C 7 new (m w ) + 0.09 C^ cw (m w ) , 



3 For consistency, the New Physics contributions entering the expression for C 7 should be taken at 
next-to-leading order in a s (mw)> i- e -j ha the radiative decay width the corresponding next-to-leading 
order New Physics matching corrections would be accounted for through C 7 rather than C 7 . 
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Table 3: Ranges of £(M) for various New Physics contributions to C-j and 
characterized hy the particles in penguin diagrams 



Class- 1 models 


e(M) 


Class-2 models 


£(M) 


neutral scalar-vectorlike quark 


i 


scalar diquark-top 


4.8-8.3 


gluino-squark (m§ < 1.37rriq ) 


-(0.13-1) 


gluino-squark (m§ > 1.37mg) 


-(1-2.9) 


techniscalar 


« -0.5 


charged Higgs-top 


-(2.4-3.8) 






left-right W-top 


w -6.7 






Higgsino-stop 


-(2.6-24) 



C 8 « -0.15 + 0.70 Cf w (m w ), 

C$ 1} » 0.48 - 2.29 C™ w {m w ) - 0.12 C 8 ncw (m w ) . (21) 

Below, we will parametrize our results in terms of the magnitude and phase of one 
of the New Physics contributions, C% ew (mw) = Kge^ or C^ ew (mw) = —K^e 117 , as well 
as the ratio 

= CT(m w ) 

where Qd — ~\- A given New Physics scenario will make predictions for these quantities 
at some large scale M. Using the renormalization group, it is then possible to evolve 
these predictions down to the scale mw- At leading order, the analogues of the relations 
© imply 

£ = £{m w ) = r£{M) - 8(1 - r) , Cf w (m w ) = r 7 C 8 ncw (M) , (23) 

where r = [a s (M)/a s (mw)] 2 ^ 3b - Here b = 11 — ^nj — 2n g is the first /5-function coefficient, 
rif = 6 is the number of light (with respect to the scale M) quark flavours, and n g — 0, 1 
denotes the number of light gluinos. For the purpose of illustration, let us consider the 
three values M = 250 GeV, 1 TeV and 2.5 TeV, which span a reasonable range of possible 
New Physics scales. We find 

f w 0.98 £(250 GeV) - 0.12 - 0.03n g 
w 0.97f(l TeV) - 0.23 - 0.03n 9 

~ 0.96 £(2.5 TeV) — 0.29 — 0.04n 9 , (24) 

i.e. £ tends to be smaller than £(M) by an amount of order —0.1 to —0.3 depending on 
how close the New Physics is to the electroweak scale. These relations will be useful for 
the discussion below. 

For simplicity, we shall restrict ourselves to cases where the parameter £ in (|2~2"D is 
real. (Otherwise there would be even more potential for CP violation.) This happens if 
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there is a single dominant New Physics contribution, such as the virtual exchange of a 
new heavy particle, contributing to both the magnetic and the chromo-magnetic dipole 
operators. Ranges of £ (M) for several illustrative New Physics scenarios are collected in 
Table |3|. They have been obtained, for simplicity, at leading order in a s and at the New 
Physics scale M characteristic of each particular model. With the help of the relations 
in (|24j), the values of £(M) can be translated into the corresponding values of £, which 
enter our theoretical expressions. Our aim here is not to carry out a detailed study of 
each model, but to give the reader an idea of the sizable variation that is possible in £. 
It is instructive to distinguish two classes of models: those with moderate (class-1) and 
those with large (class-2) values of |£|. It follows from (|2T|) that for small positive values 
of £ it is possible to have large complex contributions to C$ without affecting too much 
the magnitude and phase of C 7 , since 

Cs _ 0.70if 8 e^ -0.15 

C 7 ~ (0.09 - 0.22£)K 8 e il& ~ 0.31 ' [ ' 

This is also true for small negative values of £, albeit over a smaller region of parameter 
space. New Physics scenarios that have this property belong to class-1 and have been 
explored in Ref. ||20|| . They allow for large CP asymmetries resulting from the C 7 -Cg 
interference term in (0). Examples are penguin diagrams containing new neutral scalars 
and vector-like quarks with charge Qd = — |, for which £(M) = 1 and hence £ ~ 0.8, and 
super symmetric penguins containing light gluinos and squarks, for which £ is negative 
and can be tuned by adjusting the mass ratio m§/ rriq. A detailed analysis of the decays 
B — > X s 7 in the latter scenario is given in Ref. for the case of real C 7 and Cg. In the 
table, we specifically consider graphs with flavor off-diagonal left-right down-squark mass 
insertions under the assumption that the squark masses are approximately degenerate. 
The gluino and squark masses are taken to lie in the intervals 150 GeV < rrig < 2.5 TeV 
and 250 GeV < rriq < 2.5 TeV, respectively. Another example is provided by models 
with techniscalars of charge | [20, |5|, which have £(M) ~ —0.5 and hence £ ~ —0.7. 



In class-1 models, the magnitude of Cg can be made almost an order of magnitude larger 
than in the Standard Model without spoiling the theoretical prediction for the B — > X s 7 
branching ratio. 

In Figure ^|, we show contour plots for the CP asymmetry in the (Kg, 7s) plane for six 
different choices of £ between | and —1, assuming a cutoff E 7 > 1.85 GeV on the photon 
energy (corresponding to 5 = 0.3). We repeat that the results for the CP asymmetry 
depend very little on the choice of the cutoff. For each value of £, the plots cover the 
region < Kg < 2 and < -yg < n (changing the sign of 7s would only change the sign 
of the CP asymmetry). The contour lines refer to values of the asymmetry of 1%, 5%, 
10%, 15% etc. The thick dashed lines indicate contours where the branching ratio takes 
values between 1 x 10~ 4 and 4 x 10~ 4 , as indicated by the numbers inside the squares. 
For comparison, we recall that the Standard Model prediction with this choice of 5 is 
close to 3 x 10 -4 , whereas the current experimental values are around 2.5 x 10~ 4 . The 
main conclusion to be drawn from Figure |3| is that in class-1 scenarios there exists great 
potential for sizable CP asymmetries in a large region of parameter space. Any point 
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Figure 3: Contour plots for the CP asymmetry A^ S1 for various class-1 models 



to the right of the 1% contour for A b (^, sl cannot be accommodated by the Standard 
Model. On the other hand, we see that asymmetries of several tens of percent^ are 
possible in certain extensions of the Standard Model. It is remarkable that in all cases 
the regions of parameter space that yield the largest values for the CP asymmetries are 
not excluded by the experimental constraint on the CP-averaged branching ratio. This 
is because to have large CP asymmetries the cross-products CiC* in (|T2|) are required to 
have large imaginary parts, whereas the total branching ratio is sensitive to the real parts 
of these quantities. Note, in this context, that the cutoff imposed on the photon energy 
strongly reduces the size of the coefficient of the potentially dangerous term proportional 
to | C 8 1 2 in (0) and thereby helps in keeping the prediction for the branching ratio at an 
acceptably low level even for large values of Kg. 

There are also scenarios in which the parameter £ takes on larger negative or positive 
values. In such cases, it is not possible to increase the magnitude of C 8 much over its 
Standard Model value, and the only way to get large CP asymmetries from the C 7 -C 8 
or C7-C2 interference terms in (|12D is to have C7 tuned to be very small; however, this 
possibility is constrained by the fact that the total B — > X s 7 branching ratio must be 
of an acceptable magnitude. That this condition starts to become a limiting factor is 
already seen in the plots corresponding to £ = — | and —1 in Figure ^. For even larger 



4 We show contours only until values ^4cp = 50%; for such large values, the theoretical expression for 
the CP asymmetry in (Q2h would have to be extended to higher orders to get a reliable result. 
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values of |£|, the Cj-Cg interference term becomes ineffective, because the weak phase 
tends to cancel in the ratio Cg/CV in (p5| ). Then the C2-C7 interference term becomes 
the main source of CP violation; however, as discussed in Section [| it cannot lead to 
asymmetries exceeding a level of about 5% without violating the constraint that the 
B — ► X s 7 branching ratio not be too small. Models of this type belong to the class-2 
category. Some examples are listed in the right portion of Table |3] and can be summarized 
as follows. 

Models with gluino-squark loops can have large negative £ if the ratio mg/m,q is 
sufficiently large. Penguin graphs in left-right symmetric models with right-handed 
couplings of the W boson to the top and bottom quarks and internal top-mass chirality 
flip have £(M) ~ £ ~ —6.7. Charged-Higgs-top penguins in multi-Higgs models always 
have £(M) < —2 because of the charge of the top quark. In the table graphs with internal 
chirality flip are considered, with charged Higgs mass lying in the range 125 GeV < 
m H~ < 2.5 TeV (where £ increases as m H - is increased). In general multi-Higgs models 
these graphs are enhanced by a power of mf/rub relative to their counterparts with 
external chirality flip. Examples are type-3 two-Higgs-doublet models |2~5J, left-right 



symmetric models pq| , |49||-[pljl, or models with additional Higgs doublets which do 
not acquire significant vacuum expectation values. In all of these examples new CP- 
violating phases can enter the penguin graphs, unlike in type-2 two-Higgs doublet models. 
Chargino-stop penguins always lead to sizable negative values of £. For simplicity, we 
have considered loops that contain a pure charged Higgsino which flips chirality. The 
superpartners of new Higgs doublets with negligible vacuum expectation values would, 
for example, be pure Higgsinos. The physical stop and Higgsino masses are varied in the 
ranges 175 GeV ^ m^m^ < 2.5 TeV and 125 GeV < m~ h < 2.5 TeV, respectively, under 
the simplifying assumption that the stop mass matrix has equal diagonal entries, m 2 , 
and equal off-diagonal (left-right) entries, /x 2 , with magnitudes satisfying |/i| 2 < \mm t \. 
Finally, large positive values of £ arise from penguin graphs with a charge — ~ scalar 
"diquark" and anti-top quark in the loop. The range of values for £(M) quoted is 
again obtained for graphs with internal chirality flip, and scalar diquark mass in the 
range 250 GeV-2.5 TeV (where £ decreases as the scalar mass increases). In general, the 
phase structure of new penguin contributions with internal and external chirality flip 
will differ in the above examples; however, since the former tend to dominate due to 
chiral enhancement of order nip/mb, where irip is the mass of the heavy fermion in the 
loop, £ will be real to good approximation. 

For a graphical analysis of class-2 models it is convenient to choose the magnitude 
and phase of the new-physics contribution C^ ew (mw) = —Kje 117 as parameters, rather 
than K$ and 7s- The reason is that for large |£| it becomes increasingly unlikely that 
Cg ew (m\y) will be large. The resulting plots are given in Figure |j. As before, the dashed 
lines indicate the acceptable range for the B — > X s ^ branching ratio. The branching-ratio 
constraint allows larger values of Cg for positive £, which explains why larger asymmetries 
are attainable in this case. For example, for £ ~ 5, which can be obtained from scalar 
diquark-top penguins, asymmetries of 5-20% are seen to be consistent with the B — > X s ^ 
bound. On the other hand, for £ w —(2.5-5), which includes the multi-Higgs-doublet 
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Figure 5: B —> X sg branching ratio (left) and charm yield in B decays (right) as a 
function of the parameters Kg and 73. There is an overall theoretical uncertainty of 
6% on the values of n c . 
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models, CP asymmetries of only a few percent are attainable, in agreement with the 
findings of previous authors [^5|, ^6], [52|. The same is true for the left-right symmetric 
WM;op penguin, particularly if one takes into account that Kr ^ 0.2 if % B > 1 TeV. 

The New Physics scenarios explored in Figure ^| have the attractive feature of a 
possible large enhancement of the magnitude of the Wilson coefficient C%. This has 
important implications for the phenomenology of the semileptonic branching ratio and 
charm production yield in B decays, through enhanced production of charmless hadronic 
final states induced by the b — > sg flavour-changing neutral current (FCNC) transition 
21, j53|. At 0(a s ), the theoretical expression for the B — > X sg decay rate is obtained 



from obvious substitutions in (|6]) to be 

r(B^J = ^M ! (26) 

Normalizing this to the semileptonic rate, we obtain for the corresponding branching 
ratio B(B — > X sg ) « 0.96 |C 8 | 2 x B(B — > X c eu). In the first plot in Figure j|, we show 
contours for the B — > X sg branching ratio, normalized to B(B — > X c ev) = 10.5%, in 
the (i^ 8 ,7g) plane. In the Standard Model, B(B — > X sg ) « 0.2% is very small; however, 
in scenarios with | C 8 1 = 0(1) sizable values of order 10% for this branching ratio are 
possible, which simultaneously lowers the theoretical predictions for the semileptonic 
branching ratio and the charm production rate n c by a factor of [1 + B(B — > X sg )] . 
The most recent value of n c reported by the CLEO Collaboration is 1.12 ± 0.05 [ f40[ . 
Although the systematic errors in this measurement are large, the result favours values 
of B(B — > X sg ) of order 10% |5l}] . This is apparent from the second plot in Figure [| 



where we show the central theoretical prediction for n c as a function of K$ and 7§. 
(There is an overall theoretical uncertainty in the value of n c of about 6% resulting 
from the dependence on quark masses and the renormalization scale.) The theoretical 
prediction for the semileptonic branching ratio would have the same dependence on K% 
and 7s, with the normalization B SL = (12 ± 1)% fixed at K 8 = p5| . A large value 
of B(B — > X sg ) could also help in understanding the v( yields in charmless B decays 
p^ . For completeness, we note that the CLEO Collaboration has recently presented 
a preliminary upper limhf] on B{B —* X sg ) of 6.8% (90% CL) |60| |. It is therefore 



worth noting that large CP asymmetries of order 10-20% are easily attained at smaller 
B — > X sg branching ratios of a few percent, which would nevertheless represent a marked 
departure from the Standard Model prediction. 



5 Dipole operators with right-handed light quarks, 
and models without CKM unitarity 

All the models listed in Table ^ can have non-standard dipole operators involving right- 
handed light-quark fields. In fact, in the absence of horizontal symmetries which impose 

The limit is increased to 8.9% if one uses the more recent charmed baryon and charmonium yields 
presented in Refs. |^o|, Q and makes use of the relative A c versus A c yields given in Ref. JE^]. 
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special hierarchies among the model parameters there is no reason why these should be 
any less important than the operators of the standard basis. We therefore briefly discuss 
modifications to our previous analysis in their presence. Denoting by and C% the 



Wilson coefficients multiplying the new operators, the expressions dl2|) , fllq) and (f2q) 
must be modified by replacing CjC* — > CiC* + C^Cf* everywhere, taking however into 
account that Qf = 0. Note that for a single dominant New Physics contribution the 
parametrization in ( p2|) for the standard dipole operators will also be valid for the new 
operators, with £ taking the same real value. Then the only change in the prediction for 
the CP asymmetry is that in the denominator of ( |12"1) the coefficient |CV| 2 is replaced by 
|CV| 2 + | | 2 . On the other hand, there are several new contributions to the prediction 
for the total B — > branching ratio, as can be seen from fll8|) . For the purpose of 
illustration, let us assume that the New Physics contributions are the same for operators 
of different chirality, i.e. C^ new (m w ) = C™ w (m w ) for i = 7,8. The results are shown in 
Figure where we explore the same range of £ values as in Figure ^. The predictions 
for the B — > X sg branching ratio are enhanced because |C§| 2 in ( p6|) is replaced by 
| C 8 1 2 + |C^| 2 , so we only consider the range < K 8 < 1.5, which covers the same values 
of B(B — > X sg ) as before. Comparing Figures |3] and |6], we observe that although there 
is a clear dilution of the resulting CP asymmetries caused by the inclusion of opposite- 
chirality operators, there is still plenty of parameter space in which the asymmetries are 
much larger than in the Standard Model. We should also point out that, if there is more 
than one significant New Physics contribution to the dipole operators, there need not be 
any dilution since the product C$C^* could develop an imaginary part, thus providing 
an additional contribution to the CP asymmetry. 

Finally, we briefly discuss what happens in models with CKM unitarity violation. 
In terms of the quantity A s defined by v u + v c + (1 + A s )v t = 0, the result for the CP 
asymmetry in (ffh generalizes to 



T(B^X sl )-T(B^X- s7 ) 



T(B^X s <y) + T(B^X s <y) 



E-,>(1-5)E™™ 



a (mA f 40 8z r i 

= -j^r { si M(l + A S )C 2 CV1 - — [v(z) + b(z, 6)} Im[(l + e s + A S )C 2 C* 7 ] 

- ~ Im[C 8 CV*] + ^ b(z, 5) Im[(l + e s + A S )C 2 C* 8 ] j . (27) 

A s parametrizes the deviation from unitarity of the 3-generation CKM matrix, which 
could be caused, for instance, by mixing of the known down quarks with a new isosinglet 
heavy quark, or by the existence of a sequential fourth generation of quarks. In principle, 
asymmetries much larger than in the Standard Model could be attained provided that 
A s has a significant weak phase. This reflects the fact that the GIM suppression is no 
longer at work if CKM unitarity is violated. However, we will now show that in plausible 
scenarios the effect of A s on the CP asymmetry is very small. In the case of mixing with 
isosinglets, existing experimental limits pTJ on the FCNC process B — > X s £ + £~ induced 
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Figure 6: Contour plots for the CP asymmetry A b ^, sl for the same choices of the 
parameter £ as in Figure ^, but including the effect of different-chirality operators (as 
explained in the text) 



by tree-level Z exchange |3| imply A s < 0.04. The impact of non-unitarity can therefore 
be safely neglected, since new contributions to the CP asymmetry would be well below 
1%. Let us, therefore, turn to the case of a sequential fourth generation with a new 
up-type quark denoted by t' . As before, we will neglect the small quantity e s , so that 
Im[A s ] is the only source of CP violation. Then the above expression can be rewritten 
in the simpler form 



A b c -; s \5) = a 27 (5) Im 



;i + a 5 )c 2 



+ a s7 Im 



C 8 



C 7 



+ a 28 (<5) Im 



;i + A g )c 2 ci 

' c* 7 



C 7 



■ (28) 



In such a scenario, the CP asymmetry is affected not only by the non-unitarity of the 



3-generation CKM matrix with A s = Vf/vt in (27]), but also by the new contributions 
of the t' quark to the Wilson coefficients C 7 and C§ at the scale mw I n analogy with 
(p0|), we have 



C 7 (m w ) = ~-[A(x t ) + A s A{x t >) 



C 8 (m w ) = --[D(x t ) + A s D(x t , 



(29) 



where Xf = (mt> (mw) / n^w) 2 ■ 111 addition, there is a modification to the evolution 
equations (0) for the Wilson coefficients C 7 and Cg, where now the last terms (those 
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involving the coefficients hi and hi) must be multiplied by — (v c + v u )/v t = (1 + A s ). 
Taking m t i = 250 GeV for the purpose of illustration, we obtain C 7 « — 0.31 — 0.34A S and 
Cs ~ — 0.15 — 0.16A S , i.e. to a good approximation we have C-j^ ~ (l+A s )Cf g 1 . This just 
reflects the fact that the functions A(x) and D(x) are slowly varying for x ^> 1. In this 
limit, however, all dependence on A s cancels in the expression for the CP asymmetry. 
As a result, there is in general not much potential for having large CP asymmetries in 
models with a sequential fourth generation. For all realistic choices of parameters, we 
find asymmetries of less than 2%, i.e. of a similar magnitude as in the Standard Model. 

6 Conclusions 

We have presented a study of direct CP violation in the inclusive, radiative decays 
B — > X s j. From a theoretical point of view, inclusive decay rates entail the advantage 
of being calculable in QCD, so that a reliable prediction for the CP asymmetry can 
be confronted with data. From a practical point of view, it is encouraging that the 
rare radiative decays of B mesons have already been observed experimentally, and high- 
statistics measurements of the corresponding rates will be possible in the near future. 
We find that in the Standard Model the CP asymmetry in B — > X s ^ decays is strongly 
suppressed by three small parameters: a s (mj) arising from the necessity of having strong 
phases, sin 2 #c ~ 5% reflecting a CKM suppression, and (m c /mb) 2 8% resulting from a 
GIM suppression. As a result, the CP asymmetry can be safely predicted to be of order 
1% in magnitude. This conclusion will not be significantly modified by long-distance 
contributions. We have argued that the latter two suppression factors are inoperative 
in extensions of the Standard Model for which the effective Wilson coefficients C 7 and 
Cg receive additional contributions involving non-trivial weak phases. Much larger CP 
asymmetries of 0(a s ) are therefore possible in such cases. 

We have presented a model-independent analysis of New Physics scenarios in terms 
of the magnitudes and phases of the Wilson coefficients C 7 and C 8 , finding that, indeed, 
sizable CP asymmetries are predicted in large regions of parameter space. Some explicit 
realizations of models with large CP asymmetries have been illustrated. In particular, we 
have shown that asymmetries of 10-50% are possible in models which allow for a strong 
enhancement of the contribution from the chromo- magnetic dipole operator. This is, 
in fact, quite natural unless there is a symmetry that forbids new weak phases from 
entering the coefficients C 7 and C 8 . We have also shown that the predictions for the CP 
asymmetry are only moderately diluted if operators involving right-handed light-quark 
fields are included in the analysis. On the other hand, we confirm the findings of previous 
authors regarding the smallness of the CP asymmetry that is attainable in two-Higgs- 
doublet models and in left-right symmetric models. Moreover, we find very small effects 
for models in which 3-generation unitarity is violated. Quite generally, having a large 
CP asymmetry is not in conflict with the observed value for the CP-averaged B — > X s 7 
branching ratio. On the contrary, it may even help to lower the theoretical prediction for 
this quantity, and likewise for the semileptonic branching ratio and charm multiplicity in 
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B decays, thereby bringing these three observables closer to their experimental values. 

The fact that a large inclusive CP asymmetry in B — > X s ^ decays is possible in many 
generic extensions of the Standard Model, and in a large region of parameter space, offers 
the exciting possibility of looking for a signature of New Physics in these decays using 
data sets that will become available during the first period of operation of the B factories 
(if not existing data sets). A negative result of such a study would impose constraints on 
many New Physics scenarios. A large positive signal, on the other hand, would provide 
interesting clues about the nature of physics beyond the Standard Model. In particular, 
a CP asymmetry exceeding the level of 10% would be a strong hint towards enhanced 
chromo-magnetic dipole transitions caused by some new flavour physics at a high scale. 

We have restricted our analysis to the case of inclusive radiative decays since they 
entail the advantage of being very clean, in the sense that the strong-interaction phases 
relevant for direct CP violation can be reliably calculated. However, if there is New 
Physics that induces a large inclusive CP asymmetry in B — > decays, it will in- 
evitably also lead to sizable asymmetries in some related processes. In particular, since 
we found that the inclusive CP asymmetry remains almost unaffected if a cut on the 
high-energy part of the photon energy spectrum is imposed, we expect that a large asym- 
metry will persist in the exclusive mode B —>■ K*j, even though a reliable theoretical 
analysis would be much more difficult because of the necessity of calculating final-state 
rescattering phases Still, it is worthwhile searching for a large CP asymmetry in 
this channel. 

Finally, it has been shown in Ref. 0] that New Physics can lead to a large time- 



dependent CP asymmetry in exclusive B° — > K* 0r y decays through interference of mixing 
and decay. Large direct CP violation would introduce hadronic uncertainties, thus com- 
plicating the analysis of this effect. However, it is interesting to note that the two 
phenomena are in a sense complementary in that to a large extent they probe different 
New Physics contributions. We have seen that direct CP asymmetries in radiative B 
decays are primarily sensitive to modifications of the Wilson coefficients of the dipole 
operators with standard chirality. On the other hand, the presence of dipole operators 
with right-handed light-quark fields, which are of negligible strength in the Standard 
Model, is crucial for obtaining time-dependent asymmetries, since these require both the 
B° and B° to be able to decay to states with the same photon helicity. 
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